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In proton-nucleus and nucleus-nucleus collisions up to central S — U interactions, 
the -Py- dependence of J/ifi production is determined by initial state parton scattering 
and pre-resonance nuclear absorption ( "normal" J/ ifr suppression). The "anomalous" 
J/ip suppression in Pb — Pb collisions must reduce the normal Pt broadening, since it 
occurs mainly in the central part of the interaction region, where also initial state parton 
scattering and nuclear absorption are strongest. We thus expect for (Pj>) in Pb — Pb 



collisions a turn-over and decrease with increasing Et- 
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The recently announced "anomalous" J/ip suppression in Pb—Pb collisions [1] has 
provided new support for the hope that colour deconfinement can be established at 
CERN-SPS energy. The estimated energy densities for such collisions do fall into the 
deconfinement regime predicted by lattice QCD [2], so that there is a theoretical basis 
for such expectations. 

J/ip suppression was predicted to signal deconfinement [3]. However, all J/ip pro- 
duction in nuclear collisions up to central S — U interactions shows only pre-resonance 
absorption in nuclear matter ("normal J/ip suppression") [4,5] and hence no evidence 
for deconfinement. In contrast, Pb — Pb collisions are found to suffer an additional 
("anomalous") suppression, increasing with centrality [1]. So far, this behaviour cannot 
be consistently accounted for in terms of hadronic comovers [5] , in spite of a number of 
attempts [6-8]. It is compatible with an onset of deconfinement [5] [9-11]. 

Up to now, the experimental studies of anomalous J/tjj suppression have addressed 
only its centrality dependence. In this note, we want to consider how such suppression 
should depend on the J/ip transverse momentum. It is well known that the trans- 
verse momenta of secondaries from hadron-nucleus collisions quite generally show a pt- 
broadening. For secondary hadrons, this is the Cronin effect [12]; a similar behaviour is 
observed also in Drell-Yan and charmonium production. The natural basis for all such 
broadening is initial state parton scattering, and it was in fact shown some time ago 
[13-16] that this describes quite well the pr-dependence observed in J/ip production 
from p — A to central S — U collisions. 

Consider J/ip production in p— A collisions, assuming gluon fusion as the dominant 
process for the creation of a cc pair. Parametrizing the intrinsic transverse momentum 
distribution f(qr) of a gluon in a nucleon as 

i r # 



nqT)= m e * p vmr (1) 

we obtain by convolution for the transverse momentum distribution F p a(Pt) of the 
resulting J/ip 

i r p| 



with 



The quantity 



f » a{Pt) =w^^\-w^^ (2) 



[Pt)pA = (q 2 T )A + (q 2 T ) P . (3) 



5 pA = (P%) pA - (P%) pp = {q 2 T ) A - (q 2 T ) p (4) 



is thus a suitable measure for the observed nuclear broadening. 

Assume now that in the passage of the projectile proton through the nuclear target, 
successive interactions broaden the intrinsic momentum distribution of the correspond- 
ing projectile gluon which will eventually fuse with a target gluon to form a J/ip [13-15]. 
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If the process of Pt broadening during the passage is a random walk, then the relevant 
parameter of the Gaussian distribution (1) becomes 

5 pA = N?5 Q , (5) 

where N^ is the average number of collisions the projectile undergoes on its passage 
through the target up to the fusion point, and do the average broadening of the intrinsic 
gluon distribution per collision. 

In nucleus-nucleus collisions, a corresponding broadening occurs for both target and 
projectile gluon distributions; here, however, measurements at fixed tranverse hadronic 
energy E? can determine the broadening for collisions at a given centrality. Hence at 
fixed impact parameter b we have 



S AB (b) = (P^)AB(b) - (P|) PP = N? B (b) So, (6) 

with N^ B (b) denoting the average number of collisions for projectile nucleons in the 
target and vice versa, at fixed b. N^ B (b) has a maximum at small b and then decreases 
with increasing 6; for a hard sphere nuclear model, it would vanish when b = Ra + Rb- 

In Glauber theory, the quantity N^ B {b)/a can be calculated parameter- free from 
the established nuclear distributions [17]; here a denotes the cross section for the inter- 
action of the nucleon on its passage through the target. We shall determine abo from 
data, so that a never enters explicitly. Once ado is fixed, the broadening by initial state 
parton scattering is given for all p —A and A—B interactions. For Drell-Yan production 
(with quarks instead of gluons in the partonic interaction), this would be the observed 
effect, since the final state virtual photon does not undergo any further (strong) in- 
teractions. A produced nascent J/ip will, however, experience pre-resonance nuclear 
absorption; this suppresses J/^'s produced early along the path of the projectile, since 
they traverse more nuclear matter and hence are absorbed more than those produced 
later. As a net result, this shifts the effective production point to a later stage. In 
p —A collisions, a Drell-Yan pair will on the average be produced in the center of the 
target. In contrast, nuclear absorption shifts the average cc production point further 
down-stream. This effectively lengthens the path for initial state parton scattering and 
hence increases the resulting broadening. 

The transverse momentum behaviour of normal J/ip production in nuclear collisions 
is thus a combination of initial state parton scattering before the production of the basic 
cc state, and pre-resonance nuclear absorption afterwards; both lead to a broadening of 
a Py-distribution. A further broadening could arise from elastic random walk scattering 
of the charmonium state itself in nuclear matter; however, such an effect will be included 
here if we fit a8o to the data. 

The essential task is thus to calculate the number of collisions per cross section, 
N c /a, for p —A and A — B interactions, taking into account the effect of pre-resonance 
nuclear absorption. 



We begin with p— A collisions. The number of collisions which the projectile nucleon 
undergoes up to the cc formation point (6, z) inside the target is given by 

N c (b,z) 



a 



AT A (b,z) = Af dz'p(b,z'), (7) 

J — oo 



with b denoting the impact parameter. The total effective number of collisions before 
cc formation is obtained by averaging over z with the nuclear density p(b, z) and the 
survival probability S A (b, z) under nuclear absorption as weights, 

' Sr[h) I dz p{b,z) Nc ^ z ^ S A (b,z)/ [ dz p{b,z) S A {b,z). (8) 



(9) 



The survival probability in Eq. (8) is given by 



S A (b, z) = exp I -{A - 1) a abs / dz'p A (b, z') \ 



where a a bs is the absorption cross section for pre-resonance charmonium in nuclear 

matter; in a recent analysis [5], it was found to be 7.3 ± 0.6 mb. Eqs. (8) and (9) lead 

to 

N c (b) AT A (b) A 



a l-exp{-(A- l)a abs T A (b)} (A-l)a abs ' 

From Eq. (10) we get 

N c (b) 



(10) 



a 

and 



N c (b) 



= Wa(P), (11) 

^"063=0 ^ 

= AT A (b); (12) 

for no pre-resonance absorption, the formation thus occurs in the center of the nucleus, 
for infinite absorption on the far surface, as expected. 

Averaging the number of collisions over the impact parameter b then gives us the 
required N c /a for p —A interactions, 

^ = J d% ^) [i- P Q (b)]l j d 2 b [1- P (6)]. (13) 

Here Po(b) = [1 — aT A (b)] A is the probability for no interaction of the projectile; for a 
hard sphere nuclear distribution, it would become 0(6 — R A ). 

The formulation of the corresponding expressions for A — B interactions is quite 
straight-forward. The cc formation point is now specified by the impact parameter 6, 
the positions (s, z) and (b — s, z') in the two nuclei, with s in the transverse plane and 
z, z' along the beam axis. The number of collisions up to the formation point becomes 



N c (b,s,z,z') 



a 



A dz A p A (s,z A ) + B dz B pB{b- s,z B ), (14) 

J — 00 J —00 



and the corresponding average number of collisions in the presence of pre-resonance 
absorption is for fixed impact parameter b given by 



Nc(b) /\o r *„.,„ ^ r U - <U „ .AC- ^ O ,, „ .A ^(M,*,*') 



<i 2 s / dz pa(s, z) l dz' psib — s, z') Sa(s, z) Ssib — s,z 

J-00 J-oo & 

/ d 2 s j dz pa(s, z) / dz' ps(6 — s, z) Sa(s, z) Ss(b — s, 2/). 



(15) 



From this, we can in turn obtain the corresponding value at fixed transverse energy Et 
in the usual fashion, 

^M = Jd 2 bP(E T ,b) [1-P (6)] ^1 j dHP{E T ,b) [1-P (6)], (16) 

by convolution with the -Ex — 6 correlation function P(Ex,b) [5]. -Po(fr) here denotes 
the probability for no interaction in A — S collisions, a generalisation of the p —A form 
used above. 

With Eq. (13) for p —A and Eqs. (15/16) for A—B collisions, we have the required 
Glauber results. Ideally, we would now use Eqs. (4/5) and p —A data to fix aSo; the 
broadening for A — B interactions would then be fully predicted. Unfortunately there 
are p —A data only for three values of A [18,19], and these have rather large errors. 
We shall therefore instead check if we can obtain a consistent description of all existing 
p— A [18,19] and A—B data [19,20], up to central S—U, in terms of a common ado. F° r 
(Pt)ppi NA3 data on J/ip production at 200 GeV beam momentum [18] give 1.23 ±0.05 
GeV 2 ; this value is confirmed in an NA38 analysis [19] as well as by a two-parameter fit 
to all data which we have performed. Using the NA3 proton-proton value, we obtain the 
best fit to the available p— A, O—Cu, O—U and S—U data with ado = 9.4 ±0.7; the error 
corresponds to 95% c. L, and the minimum x 2 /d-f- is 1.1. In the Table, we show the 
experimental results together with the broadening as obtained from our calculations, 
using the mentioned a a b s = 7.3 ± 0.6 mb for the pre-resonance absorption cross section. 
The behaviour of {Pj)su as function of Et is shown in Fig. 1. It is seen that initial 
state parton scattering and pre-resonance absorption indeed account quite well for the 
observed Et dependence. 

We now turn to Pb—Pb collisions; the corresponding "normal" transverse momen- 
tum behaviour is shown in Fig. 2* Its basic feature remains the monotonic increase of 
(P T ) with Et, even though the collision geometry makes this slightly weaker for Pb—Pb 
than for S — U interactions [21]. 

The onset of anomalous suppression results in a striking modification of this pat- 
tern. If the J/'i/Vs in the hot interior of the medium produced in the collision are 
suppressed, then this will reduce their contribution from the part of phase space leading 



* A similar form was recently obtained from initial state parton scattering only, 
neglecting pre-resonance absorption [21]. Since such absorption can be simulated by 
choosing a somewhat larger So, the overall pattern remains comparable. 



to the most broadening. To illustrate the effect, we assume suppression by deconfine- 
ment; in this case the result is readily calculable [5,9,16]. To be specific, we assume 
deconfinement to start once the interaction measure [5] 

_ JV C (&, s) 

«(M) = 17 

N w (b,s) 

passes a certain critical value k c ; N c and N w are the number of collisions and of wounded 
nucleons, respectively, as obtained from Glauber theory. To calculate the effect of such 
a melting in the hot center, the integrations in Eq. (15) just have to be constrained to 
the cool outer region by introduction of 0(k c — n(b, s)). In Fig. 2 we have included the 
resulting patterns for several value of k c ; the change of behaviour due to anomalous 
J/ip suppression is qualitatively evident. 

To obtain a quantitative prediction for the forthcoming Pb—Pb results, we take into 
account the mixed (40%/60%) origin of the observed J/z/>'s from intermediate x c 's and 
direct ip production. In this case, the overall anomalous suppression had been found 
[5] to be quite well reproduced with K* ~ 2.3 and nf ~ 2.9. Using these values and 
the mentioned x/V composition, we compare in Fig. 3 the anomalous PV-dependence 
of (P T ) to the normal form. Deconfinement should thus lead to an onset of anomalous 
behaviour in the Pr-dependence of J/ip production just as it occurs in the integrated 
production rates. - In contrast, the Pt broadening observed for Drell-Yan production 
in p — A and A — B collisions [18-20] should continue in its "normal" fashion also for 
Pb — Pb interactions. 

The observation of such anomalous transverse momentum behaviour for J/ip sup- 
pression would thus constitute another piece of evidence for the onset of new physics in 
Pb — Pb collisions. In particular, its occcurrence in Pb — Pb collisions, with a normal 
Pr-dependence in S — U interactions, would again rule out a hadronic explanation of 
anomalous J/tf) suppression: any absorption by hadronic comovers must be present in 
continuously varying degrees in all nucleus-nucleus collisions. On the other hand, as we 
have seen here, such anomalous P^-behaviour would arise naturally in case of deconfine- 
ment. It should disappear only for Pt values so high that gluon fragmentation replaces 
gluon fusion as dominant J/ip production mechanism or when the relative J/^-gluon 
momentum becomes too large for an effective J/tf) dissociation (for Pj- > 5 — 10 GeV) 
[22]. 



Table Caption 

Experimental and theoretical values of (P T ) for p —A and A — B interactions. Data for 
p - p and p -Pt: NA3 [18]; the other from NA38 [19,20]. 

Figure Captions 

Fig. 1: Pt broadening in S — U collisions; diamonds from [19], stars from [20]. 

Fig. 2: Normal and anomalous Pt behaviour for Pb — Pb collisions; the anomalous 
behaviour is calculated for deconfinement at the given critical values of the interaction 
measure k c . 

Fig. 3: The Pt behaviour for deconfinement in Pb — Pb collisions for a 40%/60% x/tf 
origin of the observed J/ip's, compared to normal behaviour. 
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< 13 
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